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Abstract 

Results obtained with the a new Monte-Carlo code for the channelling of 
855 MeV electrons along the crystallographic plane (110) in a bent Silicon 
crystal are presented. The definitions of the dechannelling length and the 
asymptotic acceptance of the channel are given in a form that is suitable 
for the Monte-Carlo procedure. The dependence of these quantities on the 
crystal bending is studied. 
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1. Introduction 

In this article we study the channelling of ultrarelativistic electrons in a 
bent planar channel using a new Monte-Carlo code. 

Channelling takes place if charged particles enter a single crystal at small 
angle with respect to a major crystallographic direction: a plane or an axes 
[l[. The particles get confined by the electrostatic potential of the crystal 
and move preferably along the corresponding crystallographic planes or axes 
following their shape. This phenomenon has been intensively studied since 
the Sixties of the last century. Recently, practical applications of channelling 
have been inducing growing interest in the scientific community. 
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In particular, the crystals with bent crystallographic planes are used 
to steer high-energy charged particle beams replacing huge dipole magnets. 
Since its appearance [2!] and first experimental verification jsl this idea has 
been attracting a lot of interest worldwide. Bent crystals have been routinely 
used for beam extraction and deflection in the Institute for High Energy 
Physics, Russia A series of experiments on the bent crystal deflection 
of proton and heavy ion beams were performed j^, [g], 0, S, [oj. The bent 
crystal method has been proposed to extract particles from the beam halo 
at CERN's Large Hadron CoUider [lo|. The possibility of deflecting positron 
m and electron [o]. Hi], 13 1 beams has been studied as well. 

Another very promising application of the channelling phenomenon is the 
crystalline undulator, — a novel source of hard electromagnetic radiation. A 
single crystal with periodically bent crystallographic planes can force chan- 
nelling particles to move along nearly sinusoidal trajectories and radiate in 
the hard x- and gamma-ray frequency range. The feasibility of such a device 
was demonstrated theoretically a decade ago fl3. [l5| (further developments 



as well as historical references are reviewed in 16|). The advantage of the 
crystalline undulator is the extremely strong electrostatic fields inside a crys- 
tal, which are able to steer the particles much more effectively than even the 
most advanced superconductive magnets. This fact allows one to make the 
period of a crystalline undulator very short (up to several microns) and, there- 
fore, the wavelength of the produced radiation can reach the (sub)picometer 
range, where conventional sources with comparable intensity are unavailable 



_17|. 

Initially, it was proposed to use positron beams in the crystalline un- 
dulator. Positively charged particles are repelled by the crystal nuclei and, 
therefore, they move between the crystal planes, where there are no atomic 
nuclei and the electron density is less then average. This reduces the probabil- 
ity of random collisions with the crystal constituents. Hence, the transverse 
momentum of the particle increases slowly and the particle travels a longer 
distance in the channelling regime. 

More recently, an electron based crystalline undulator was proposed [18|, 



19| . On one hand, electrons are less preferable than positrons. Due to their 
negative charge, the electrons are attracted by the lattice ions and, therefore, 
are forced to oscillate around the crystal plane in the process of channelling. 
The probability of collisions with crystal constituents is enhanced. Thus, the 
dechanneling length is smaller by about two orders of magnitude in compar- 
ison to that of positrons at the same conditions. On the other hand, the 
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electron beams are easier available and are usually of higher intensity and 
quality. Therefore, from the practical point of view, electron based crystalline 
undulator has its own advantages and deserves a thorough investigation. 

To develop a comprehensive theory of the crystalline undulator, one has to 
understand the behaviour of the electron in a bent crystal channel. To study 
this problem, we have developed a new Monte-Carlo code that allows us to 
simulate the particle channelling and analyse their trajectories. In contrast to 
other channelling codes 20|, l2l|, |22| , our algorithm does not use the continuous 
potential approximation. This novel feature is especially beneficial in the case 
of negatively charged projectiles, which channel in the vicinity of the atomic 
nuclei, where the continuous potential approximation becomes less accurate. 

In the present paper, we restrict our analysis to channels with constant 
curvature. Our results will help one to estimate the reasonable parameters 
of periodically bent crystal channels and, therefore, will facilitate future sim- 
ulations and experimental study of the electron-based crystalline undulator. 
Additionally, the results can be useful to developers of crystalline devices for 
extraction, bending and coUimation of electron beams. 

We study the channelling of 855 MeV electrons in a straight and bent 
single crystal of Silicon along the plane (110). This plane is chosen for the 



analysis because it can be deformed by growing of Sii-ajGca; crystals [23 
with a periodically varying Ge content x [ij, 25|. The electron energy corre- 
sponds to the conditions of the channelling experiments at Mainz Microtron 
(Germany) (26| . 



2. The Physical Model and the Simulation Algorithm 

In this section we review the basic ideas of the physical model and the 
algorithm implemented in our Monte Carlo code. A more detailed description 
can be found in f27| . 

The model is intended for studying the interaction of ultrarelativistic 
projectiles with single crystals. Due to the high speed of the projectile, 
its interaction time with a crystal atom is very short. The motion of the 
atomic electrons as well as the thermal motion of the atomic nuclei during 
the interaction time can be neglected. Therefore, the crystal is represented as 
a collection of static charges. The atomic nuclei 'frozen' at random positions 
in the vicinity of nodes of the crystal lattice. The probability distribution 
of the position of the nucleus relative to the node is approximated by a 
three dimensional normal distribution with the variance equal to the squared 
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amplitude of thermal vibrations of the crystal atoms a{T). We used the 
value a(T) = 0.075 A in our calculations, which corresponds to the room 



temperature [28[. The probability distributions of the electrons, given by 
squared absolute value of the wave function of the atom. Instead of the 
exact wavefunction, we use a spherical symmetric distribution that on average 
reproduces the electrostatic potential of the atom in Moliere's approximation 

m. 



The code performs 3D simulation of the projectile motion in the crystal. 
The coordinate frame is defined in such a way that the axis z coincides with 
the beam direction. The plane (110) of the crystal is parallel to the coordinate 
plane (xz) and, consequently, it is orthogonal to the axis y. The interaction 
of the projectile with a crystal constituents is considered as a classical scat- 
tering in a Coulomb field of a static point-like charge. The trajectory of 
the projectile is modelled using the following algorithm. The particle moves 
along straight line segments between the points where its projection on the 
axis z coincides with the z coordinate of a crystal constituent: an electron 
or a nucleus. At this point the transverse momentum of the projectile is 
changed by the quantity 

Apl = -gpgt— • (1) 

Here qp and qt are electric charges of the projectile and the target constituent, 
respectively, f± is the vector connecting the projectile and the target con- 
stituent (it is parallel to the axis (xy)), v ^ c is the projectile speed. Then 
the projectile is propagated further along a straight line segment until its 
z coordinates coincides with that of the next crystal constituent and a new 
modification of the momentum is performed. The ionization losses of high 



energy electrons or positrons in matter are very small [30|. The radiation 
losses become noticeable starting from the projectile energies of several tens 
GeV. In our calculations they are also neglected. Therefore, the absolute 
value of the projectile momentum remains unchanged during the simulation. 

A crystal constituent is taken into account if it belongs to the lattice node 
located within a cylinder of the radius 40aTF around the particle (otf is the 
Thomas- Fermi radius of the atom). Initially, the axis of the cylinder is the 
straight line along the direction of the projectile momentum at the point of 
entering the crystal. The length of the cylinder Azc is approximately 200 A. 
When the particle approaches the end of the cylinder, a new cylinder is built 
as an extension of the old one but along the direction of the new particle 
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momentum. The procedure continues until the end of the crystal is reached. 
As result, the cylinders form a 'pipe' filled by the crystal lattice and the 
particle channels inside it. 

In the case of a bent crystal with channel shape defined by the function 
yB{z), the above procedure is modified in the following way. If (xp, Up, Zp) are 
the coordinates of the projectile, and (xc. He, ^c) is the endpoint of the cylinder 
axis: = Zp + Azc, Xc = Xp + Az^Px/Pz, and Vc = Up + Az^Py/p^. First, an 
undistorted lattice is built around the straight line segment connecting the 
points (xp. Hp — yB{zp), Zp) and (xc, i/c — 1/b(-2c), -^c)- Then each lattice node is 
shifted according to the bending profile: (xn,yn,-2n) — ^ {xn,yn + l/B(-2n), -^n)- 
As a result, the trajectory becomes surrounded by a 'pipe' filled by the bent 
crystal lattice. 

It is convenient to characterise the curvature of the bent crystal by the 
centrifugal parameter C. Let U{v) be the potential energy of the projectile 
in the field of crystal constituents averaged over the surface y = v + yB{z) . 
Then the centrifugal parameter is defined as 

C = (2) 

^ max 

Here Fd. is the centrifugal force acting on the projectile in the bent channel 
and ?74iax is the maximum derivative of the average potential energy with 
respect to t;, i.e. the maximum force acting on the projectile in the average 
potential. The value C = corresponds to the straight crystal and C = 1 is 
the critical value at which the potential barriers between the crystal channels 
disappears. 

The centrifugal force is related to the projectile energy E and to the bend- 
ing radius Rc via Fd. = E/Rq. Therefore, expression ([2]) can be rewritten 

as 

^ = Rlj^- 

In the case of a channel with a constant bending radius, the shape is 
defined by the equation 



y^iz) = Rc- jRh - z\ (4) 



The channel boundaries for a straight crystal are defined as the planes 

where yfjp 



satisfying the equation y = y^ , where y^^l is the coordinate of the potential 
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energy maximum between the n-th and {n + l)-st channel. For C ^ these 
planes become bent surfaces following the shape of the channel: 

y = yiol + yB{z). (5) 

In this case, the channel boundaries do not coincide with the maxima of the 
potential energy which is modified by centrifugal effects (see figured]). 

3. Simulations 

The calculations were performed for E=855 MeV electrons channelling 
along the plane (110) for twelve different values of C in the range from to 
1. 

Initially, the projectiles had zero transverse momentum. This corresponds 
to the ideal case of a zero-emittance beam entering the crystal strictly par- 
allel to the coordinate axis z. The transverse position of the projectile at 
the entrance of the crystal was chosen randomly, homogeneously distributed 
along the channel width. Then the trajectory of the particle was simulated. 
The simulation of the trajectory was terminated if the particle went through 
the crystal, z > L^r {Lcr is the crystal dimensions along the beam direc- 
tion), or if the deviation of the projectile from its initial direction became 
too large: Ip^l/pz > IOO/7 (here p± and p are respectively the transverse 
and the longitudinal momenta and 7 is the Lorentz factor of the projectile.) 

Then the simulated trajectories were analysed and the segments corre- 
sponding to channelling and dechannelling mode were determined. At the 
entrance point of the crystal, the projectile was assumed to be in the chan- 
nelling mode. Regardless of the value of its transverse energy, the particle 
was considered to be channelling until it crossed the channel boundary 

The crystal dimensions along the beam direction was Lcr = 120 fim for 
a straight channel. For nonzero C, L^r was chosen so that very few or no 
particle remain in the channelling mode at the exit from the crystal. The 
number of the simulated trajectories was 100000 for the straight channel and 
25000 for each nonzero value of C. 

4. Definition of the Dechannelling Length and the Asymptotic Ac- 
ceptance of the Channel 

To make a quantitative assessment of the particle dechannelling process, 
one needs a definition of the dechannelling length that would be suitable for 
the Monte Carlo approach. 
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Figure 1: The average potential energy of the projectile in the field of the crystal planes 
modified by the centrifugal effects. The centrifugal parameter C is defined by equation 
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Let Zdi be the point of the first the channel boundary by the projectile. 
We define the quantity Ncho{z) as the number of trajectories for which 2:^1 > 
z, i.e. this is the number of particles that passed the distance from the 
crystal entrance to the point z in the channelling regime and dechannel at 
some further point. The length L{z) is the average distance from the point 
z to the first dechannelling point: 

The sum in the numerator of is taken over those projectiles for which z^i > -2- 
Generally speaking, L{z) depends not only on z, but also on the angular 
distribution of the particles at the crystal entrance. Nonetheless, as it will be 
shown below, the kinetic theory of channelling suggests that, at sufficiently 
large z, Li^z^ reaches the asymptotic value that depends neither on z nor on 
the initial angular distribution. 

From the solution of the diffusion equation (see, for instance, formula 



(1.38) in |31i), one can obtain the following expression for iVcho(-z) 

00 

A^cho(2;) = A^o J] exp (-2;/L, ) . (7) 
i=i 

Here only coefficients Aj depend on the initial angular distribution of the 
particles, while the lengths depend exclusively on the properties of the 
crystal channel and the energy, charge and mass of the projectile. 

The dechannelling length is defined as the largest of the parameters 
Lj in d?]). The corresponding term dominates the asymptotic behaviour at 

Ncw{z) X No Ad exp {-z/La) . (8) 

The expression (jH]) for L{z) has the following counterpart in the kinetic 
theory 



00 



m = I dz,,{z,, - (9) 

Substituting ([8]) into ([9]) demonstrates that, indeed, the coefficient cancels 
out and L(z) becomes equal to in the asymptotic region. 

Although the diffusion equation was solved in [31] for positively charged 
projectile in harmonic potential approximation, the exponential asymptotic 
behaviour of Ncho{z), and, consequently, a constant asymptotic value of L{z) 
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is a more general result. As it will be shown in the next section, our simula- 
tions demonstrate that it is also valid for electrons. 

Hence, in our Monte Carlo procedure the dechannelling length is de- 
fined as the asymptotic value of L{z) (Q in the region where it ceases to 
depend on z. 

Due to the exponential asymptote (E]), the quantity 

A(z) = %^exp(./L,) (10) 

should approach a constant value at sufficiently large z. We shall call 
the quantity the asymptotic acceptance of the channel. The asymptotic 
acceptance depends on the initial transverse momentum distribution of the 
particles in the beam. 

5. Analysis of the Results 

The ratio Ncho{z) / Nq of the particle that never crossed the channel bound- 
ary is shown in Figure [2] as function of z for several values of C. This fraction 

1 ' ' ' ' ' ' ' =1 



0.8 r 




z (jum) 



Figure 2: Fraction Ncho{z)/No of the particles that stay in the same channel from their 
entrance into the crystal for different values of the centrifugal parameter C as function of 
the penetration depth z. The dashed lines show the corresponding exponential asymptotes 
(X exp(— z/Ld). 
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decreases rather fast and, as it was expected, has an exponential asymptotic 
behaviour. 

The quantity L{z) for the same values of the centrifugal parameter 
are plotted in figure |3l Indeed, L{z) becomes constant (within the statistical 




Figure 3: The quantity L{z) ([5]) that becomes equal to the dechannelhng length at large 
penetration depth z. The thin lines show the statistical errors. The straight dashed lines 
correspond to the asymptotic values. 

errors) at large z corresponding to exponential behaviour of the curves of 
figure [21 The asymptotic value for the straight channel is Ld{C = 0) = 
8.464 ± 0.057 /xm. 

The dependence of the ratio L^i^C) / L^{C = 0) on the centrifugal param- 
eter is shown in figure HI 

Surprisingly, the dechannelhng length remains nonzero even at C = 1, 
where the potential barriers between the channels disappear and therefore 
the transverse motion of the particle is no more constrained by the interplanar 
potential. The reason of such behaviour is likely the random scattering of 
the particles by the crystal constituents. There is a nonzero probability 
that the projectiles is scattered in such a way that it trajectory remains 
within the channel boundaries ([5]). The fraction of such particles decreases 
exponentially. Therefore, the curve for C = 1 in figure [3] behaves in a similar 
way as the for C < 1 having a constant asymptote, which corresponds to a 
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Figure 4: The ratio Ld{C) / Ld{C = 0) of the dechanneUing length in bent channels to 
the one in the straight channels as function of the centrifugal parameter C. At C < 0.3 
the ratio can be parametrised by the function (1 — C)^'^ (shown by the solid line). For 
comparison, the behaviour (1 — C)^ obtained in the diffusion theory in the parabolic 
potential approximation is shown by the dashed line. 
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nonzero dechannelling length. 

It is seen from figure SJ that the ratio L^iC) / L^iC = 0) follows the law 
(1 — C)^'^ for C < 0.3. At C > 0.3, it decreases less steeply. This can be 
interpreted as the effect of the random scattering that becomes comparable 
to the guiding effect of the interplanar potential at C ~ 0.3. 

The quantity A{z) defined by ffTU]) is plotted in figure O This quantity 
indeed approaches a constant value in the same region of z where the quantity 
L{z) does. This value is by definition the asymptotic acceptance of the 
channel for an ideally parallel beam. For a straight channel, the value Ad{C = 
0) = 0.95 ± 0.01 is obtained. 
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Figure 5: The quantity A{z) (fTO|) that becomes equal to the asymptotic acceptance of 
the channel at large penetration depth z. The thin lines show the statistical errors. The 
straight dashed lines correspond to the asymptotic values A^- 

The dependence of the ratio A^^i^C) / A^^C = 0) on the centrifugal param- 
eter is shown in figure [6l It is seen that the behaviour of the ratio is changed 
at C = 0.3: it falls as (1 — C)^ at 0.1 < C < 0.3 then becomes nearly con- 
stant at 0.3 < C < 0.9. Then it jumps up at C = 1. At very small C, the 
behaviour is close to 1 — C law, which coincides with the one obtained in the 
diffusion approach for a parabolic potential. 
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Figure 6: The ratio Ad{C)/Aii{C = 0) of the asymptotic acceptance of bent channels 
to the one in the straight channel as function of the centrifugal parameter C. At 0.1 < 
C < 0.3 the ratio can be parametrised by the function (1 — C)^ (shown by the solid line). 
For comparison, the behaviour (1 — C) obtained in the diffusion theory in the parabolic 
potential approximation is shown by the dashed line. 
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6. Conclusion and Discussion 

We presented the Monte-Carlo study of electron channelling in the (110) 
channel of a bent Sihcon crystal. The calculation were done for the beam 
energy of 855 MeV for several value of the centrifugal parameter. 

According to our simulation, the number of channelling electrons rapidly 
decreases with the penetration depth z and quickly approaches an exponen- 
tial asymptote. Similar behaviour was previously seen in the kinetic theory 
of channelling in the case of positively charged projectile. 

Wc formulated definitions of the dechannelling length and the asymp- 
totic acceptance of the channel that are suitable for application within the 
Monte Carlo approach. Our definition is consistent with the one previously 
used in the kinetic theory of channelling. 

We calculated the dechannelling length for the set of values of the cen- 
trifugal parameter. The dechannelling length decreases as (1 — C)^'^ at 
< C < 0.3. At larger C the guiding effect of the channel becomes so 
weak that the random diffusion of the particles along the bent channel is 
comparable to channelling effect. This results into the change of the be- 
haviour of L(j(C). Similarly, the behaviour of the asymptotic acceptance 
A^{C) changes at C = 0.3 from (1 — C)^ decreasing to nearly constant and 
then increasing behaviour. 

The obtained result can be used for choosing the optimum bending of the 
crystals used for deffection and collimation of electron beams. 
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